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Abstract 

The direct searches for Superymmetry at colliders can be complemented 
by direct searches for dark matter (DM) in underground experiments, 
if one assumes the Lightest Supersymmetric Particle (LSP) provides 
the dark matter of the universe. It will be shown that within the Con- 
strained minimal Supersymmetric Model (CMSSM) the direct searches 
for DM are complementary to direct LHC searches for SUSY and Biggs 
particles using analytical formulae. A combined excluded region from 
LBC, WMAP and XENON 100 will be provided, showing that within the 
CMSSM gluinos below 1 TeV and LSP masses below 160 GeV are ex- 
cluded (mi/2 > ^OOGeV) independent of the squark masses. 

1 Introduction 

Supersymmetry (SUSY) is the leading candidate for physics beyond the SM, since the Lightest Super- 
symmetric Particle (LSP) has all the properties expected for the Weakly Interacting Massive Particles 
(WIMPS) of the dark matter (THE], which is known to make up more than 80% of the matter in 
the universe [4J. Unfortunately, no supersymmetric particles have been observed so far, even at the 
highest energy of the LHC [SHE], nor have WIMPS been observed beyond doubt in elastic scat- 
tering of WIMPS on nuclei, as pursued with direct searches in underground detectors, like CDMS, 
EDELWEISS and XENON100, which give upper limits on the elastic scattering cross section typically 
between 10 -43 and 10 -44 cm 2 [HI [19]. The direct DM detection experiments, the search for SUSY 
and the search for Higgs particles all determine complementary excluded regions in the parameter 
space of supersymmetric models, like the Constrained Minimal Supersymmetric Model (CMSSM), see 
[201 El] and reviews, e.g. [22H25] . In this model all the arguments in favour of Supersymmetry, like the 
unification of the coupling constants at the GUT (Grand Unified Theory) scale with SUSY masses in 
the TeV range [26] and electroweak symmetry breaking (EWSB) by radiative corrections [23 EB] , are 
implemented. In the CMSSM one furthermore assumes that the masses of spin (spin 1/2) particles 
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are unified at the GUT scale with values mo(m 1 / 2 )- So the many parameters of SUSY models are 
reduced to only 4: the two mass parameters mo, m^ 2 and two parameters related to the Higgs sector: 
the trilinear coupling at the GUT scale Ao, and tan/3, the ratio of the vacuum expectation values 
of the two neutral components of the two Higgs doublets. Electroweak symmetry breaking (EWSB) 
fixes the scale of fi, so only its sign is a free parameter. The positive sign is taken, as suggested by 
the small deviation of the SM prediction from the muon anomalous moment, see e.g. |29| . 

Within the CMSSM the direct searches for SUSY particles (" sparticles" ) at the LHC, the direct 
DM searches and the relic density, as obtained from cosmological observations, are related and one 
can combine them to see which region of the supersymmetric parameter space is excluded, if one 
includes all constraints. Such combinations have been pursued by many different groups either using a 
frequentist approach by maximizing a likelihood or using random sampling techniques of the parameter 
space, see e.g. [301 - HT] and references therein. 

The sampling techniques are dependent on the prior, which leads to an additional, non-quantifiable 
uncertainty in the excluded or allowed regions, see e.g. [42J for a recent discussion and references 
therein. We believe this uncertainty is due to the high correlations between three of the four pa- 
rameters, as we discussed in two previous papers: mo is highly correlated with tan (3 because of the 
relic density constraint, which requires large tan/3 in most of the parameter space except for the 
narrow co-annihilation regions at low and large mo [E] and the trilinear coupling is highly correlated 
with tan/3 because of heavy flavour constraints, notably the B® —> constraint [H]. Note that 

Bg — » /i + n~ is not only proportional to tan 6 /3, so it tends to become above the present upper limit 
for large tan/3, but it can be strongly reduced by the appropriate value of the trilinear coupling to 
values even below the SM value by negative interferences [44J. 

Such strong correlations lead to likelihood spikes in the parameter region, where three of the four 
parameters have to have specific correlated values. Although the likelihood of such narrow regions is 
high, they are either not found in methods based on stepping techniques or their probability is given 
a different weight because of its "low posterior mass". To cope with the strong correlations we use a 
multistep fitting technique, defined by fitting the parameters with the strongest correlation first, i.e. 
we fit first tan /3 an Aq for each pair of the mass parameters mo and m^ 2 by minimizing the x 2 with the 
program Minuit[45j or by using a Markov Chain Monte Carlo or multinest fitting technique to find the 
maximum likelihood, while integrating over the nuisance parameters. Both the x 2 minimization and 
Markov Chain sampling of the parameter space give practical identical results in this multistep fitting 
technique and no prior dependence has been found, since for each point in the (mo,my 2 ) grid there is a 
unique solution for tan /3, mainly from the relic density constraint, and Aq can be fine-tuned to further 
minimize x 2 - Furthermore the multistep fitting technique is fast, since initially only two parameters 
are fitted for each point in the (mo,mi^) grid. Once the SUSY parameters have been fitted, one 
should also vary the SM parameters or marginalize over them, like the top and bottom mass and the 
strong coupling constant. However, these are highly correlated with the SUSY parameters, so if one 
repeats the fit with different values of the SM parameters or marginalizes over them one usually finds 
the same fit probability with slightly different values of the SUSY parameters. The SM parameters 
are given in the Particle Data Book [46J; we use m^ e = 172.5 ± 1.3 and m h {m h ) MS = 4.25 ± 0.2 GeV 
for the heavy quark masses and a s = 0.1172 zb 0.02 for the strong coupling constant. 

All observables discussed below were calculated with the public code micrOMEGAs 2.4.1 [17} H8] 
combined with Suspect 2.41 as mass spectrum calculator [49J. 

The purpose of this letter is twofold: we show that the data from the LHC, WMAP and XENON100 
lead to complementary excluded regions and discuss the reasons why by giving analytical approxima- 
tions. We furthermore combine the newest data and show that a lower limit of m^ 2 of 400 GeV can 
be obtained independent of mo, which implies a lower limit of 160 GeV on the LSP and 1 TeV on the 
gluino. Our results differ from similar analysis referred to above and we discuss possible reasons. We 
start by discussing the observations and the excluded regions of each observation separately. 
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Figure 1: Examples of diagrams for SUSY particle production via strong interactions (top rows for 
gg, qq and gq, respectively) and electroweak interactions (lowest row). 



2 Excluded region by direct searches for SUSY at the LHC 

In proton-proton collisions strongly interacting supersymmetric particles can be produced by the 
main diagrams shown in the first three rows of Fig. [TJ while the main diagrams for the electroweak 
production are shown in the last row. The corresponding cross sections are shown in Fig. [2] for a 
centre-of-mass energy of 7 TeV. One observes that the cross section for the "strong" production of 
qq and gq ) are large for low values of mo and m^ 2 , the gluino production gg is strongest at small 
values of m^ 2 and the electroweak production of gauginos starts to increase at large values of mo- 
The reason for the increase of the electroweak production at large mo is the decrease of the Higgs 
mixing parameter /x, as determined from EWSB, which leads to a stronger mixing of the Higgsino 
component in the gauginos and so the coupling to the weak gauge bosons and Higgs bosons increases, 
thus increasing the amplitudes for the diagrams in the last row of Fig. [T] 

The physics is simple: the Higgs masses at the GUT scale start with a value of yj ttiq 2 + ji 2 
and EWSB requires at least one of them to become negative at the electroweak scale by radiative 
corrections, mainly from the Yukawa couplings of the third generation. This is only possible if the 
starting value is not too high, so a large value of mo has to be compensated by a low value of /i. This 
is demonstrated in Fig. [3j which shows that n becomes small and the Higgsino component in the 
lightest neutralino becomes large for large values of mo- 

The strong production cross sections are characterized by a large number of jets from long decay 
chains and missing energy from the escaping neutralino. These characteristics can be used to efficiently 
suppress the background. For the electroweak production, both the number of jets and the missing 
transverse energy is low, since the LSP is not boosted so strongly as in the decay of the heavier 
strongly interacting particles. Hence, the electroweak gaugino production needs leptonic decays to 
reduce the background, so these signatures need more luminosity and cannot compete at present with 
the sensitivity of the strong production of squarks and gluinos. 

The total cross section for the strongly interacting particles are shown in Fig. [4] together with 
the excluded region from direct searches at the LHC for SUSY particles. One observes that the 
excluded region (below the solid line) follows rather closely the total cross section, indicated by the 
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Figure 2: Cross sections for SUSY particle production for the diagrams shown in Fig. [l| clockwise 
via strong interactions (gg, gq and qq, respectively) and electroweak interactions. 



colour shading. From the colour coding one observes that the excluded region corresponds to a cross 
section limit of about 0.1 -0.2 pb. The excluded region was obtained by combining the ATLAS [5j 
and CMS [6] limits in the following way. Since the excluded region follows rather closely the total 
cross section the x 2 contribution has been parametrized as a^ ot /a^p where can be determined 
by requiring that for each point in the mo, m^ 2 plane the x 2 value corresponds to an UP-parameter 
°f X 2 ~ Xmin ~ 5.99, which corresponds to a 95% confidence level for a two-dimensional parameter 
spac^j [45J. Determining the parameter cr 2 ff for each point of the contour and each experiment 
independently implies that our contours are identical to the contours given by the experiments. Then 
the x 2 contributions are simply added, assuming no correlation between the experiments. If the LHC 
data is combined with cosmological and electroweak data, the fitted values of tan /3 and the trilinear 
coupling vary, while the LHC limits are given for fixed values tan /3 = 10 and Aq = 0. In our definition 
of x 2 the cross section variations as function of these parameters are correctly taken into account under 
the assumption that the efficiency does not vary with these parameters, which is a good approximation 
for the hadronic searches. Therefore we only consider limits from LHC data based on jets and missing 
energy and do not include the less sensitive limits from leptonic data. 

The drop of the excluded region at large values of mo is due to the fact that in this region the 
squarks become heavy, which means that the contributions from the diagrams in the second and 
third row of Fig. [T] start to diminish. Here only higher energies will help and doubling the LHC 
energy from 7 to 14 TeV, as planned in the coming years, quickly increases the cross section for gluino 

2 This UP-parameter is for a two-sided confidence interval. However, several observables, like the LHC limits, refer 
to a single-sided lower limit. For a single-sided limit the UP-parameter would be 4.61. We keep conservatively an 
UP-parameter of 5.99, which for a single-sided limit would correspond to a 97.5% confidence level Because of the steeply 
falling cross section the difference between 95 and 97.5% is small for the exclusion line in the mo, rai/2 plane. 
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Figure 3: Values of fi and 7Vi3 in the mc^m^ plane. The strong decrease of \i and corresponding 
increase of 7Vi3 at large values of mo is largely determined by the EWSB constraint. Note 
that the right bottom corner is typically excluded by EWSB, if the trilinear coupling Aq 
and tan /3 are fixed. However, if all parameters are left free in the fit, EWSB can be fulfilled 
in the whole region shown. 



production by orders of magnitude, as shown in the right panel of Fig. [4| The expected sensitivity 
at 14 TeV, plotted as an exclusion contour in case nothing will be found, assumes the same efficiency 
and luminosity (slightly above one fb _1 per experiment) as at 7 TeV. 

These limits can be translated to squark and gluino masses as follows. The squarks have a starting 
value at the GUT scale equal to mo, but have important contributions from gluinos in the colour field, 
so the squark masses are given by m~ « + 6.Qmi/2 2 , as was determined from the renormalization 
group equations [50] . Similarly the gluino mass is given by 2.7m 1 / 2 - The term proportional to 
in the squark mass corresponds to the self-energy diagrams, which imply that if the "gluino-cloud" 
is heavy, the squarks cannot be light. This leads to the regions indicated as not allowed in Fig. [5j 
Note that these regions are forbidden in any model with a coupling between squarks and gluinos, so 
they are not specific to the CMSSM. Squark masses below 1.1 TeV and gluino masses below 0.62 TeV 
are excluded for the LHC data at 7 TeV, as shown in the left panel of Fig. [5j Expected sensitivities 
for higher integrated luminosities at 7 and 14 TeV have been indicated as well. One observes that 
increasing the energy is much more effective than increasing the luminosity. At 14 TeV squark masses 
of 1.7 TeV and gluino masses of 1.02 TeV are within reach with one fb _1 per experiment, as shown in 
the right panel of Fig. [5| 



3 Excluded region by the relic density 

The observed relic density of DM corresponds to fi/i 2 = 0.113 ± 0.004 [4J, which is about a factor six 
higher than the baryonic density. This number is inversely proportional to the annihilation cross sec- 
tion. The dominant annihilation contribution comes from A-boson exchange in most of the parameter 
space, if one excludes the narrow co-annihilation regions [33]. The cross section for % + % A bb 
can be written as: 

M^rnl tan 2 /3 (N 31 sin /3 - 7V 4 i cos /3) 2 (N 2 i cos 9 W - N u sin 9 W ) 2 

(4M 2 -M 2 ) 2 + M 2 F 2 l ' 

where the elements of the mixing matrix in the neutralino sector define the content of the lightest 
neutralino 

|X?> = Nu\B ) + N 12 \W$) + N^Hi) + N U \H 2 ). (2) 



sin 4 26 w M§ 
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Figure 4: Left: Total production cross section of strongly interacting particles in comparison with 
the LHC excluded limits for 7 TeV. Here the data from ATLAS and CMS were combined 
and correspond to an integrated luminosity of 1.3 and 1.1 fb _1 , respectively. One observes 
that a cross section of 0.1 to 0.2 pb is excluded at 95% confidence level. Right: the cross 
sections at 14 TeV and expected exclusion for the same limit on the cross section as at 7 
TeV (left panel). 




Figure 5: As in Fig. |4| but the excluded region is translated into the m^rrig plane. The red 
area corresponds to excluded regions for an integrated luminosity slightly above 1/fb; the 
expectations for higher luminosities have been indicated as well. Note that the regions 
indicated as not allowed are not a feature of the CMSSM, but are valid in any model, 
which has a coupling between squarks and gluinos, which lead to self-energy diagrams of 
the squarks. These imply that if the " gluino-cloud" is heavy, the squarks cannot be light. 

The correct relic density requires < av >= 2 • 10 -26 cm 3 /s, which implies that the annihilation cross 
section a is of the order of a few pb. Such a high cross section can be obtained only close to the 
resonance. Actually on the resonance the cross section is too high, so one needs to be in the tail of 
the resonance, i.e. vtla ~ 2.2m x or vtla ~ 1.8m x . So one expects vtla oc m^ 2 from the relic density 
constraint. This constraint can be fulfilled with tan/3 values around 50 in the whole mo,m^ 2 plane, 
except for the narrow co-annihilation regions, as we showed previously [43J. The results were extended 
to larger values of mo, as shown in the left panel of Fig. [6j The rather strong limits on the pseudo- 
scalar Higgs boson from LHC [5lJ[52], especially at large values of tan/3, lead then to constraints on 
mi/2 of about 400 GeV for intermediate values of mo, as shown in the right panel of Fig. [6j Note that 
the top left corner (white) is not allowed, because in this region of small mo and large m^ 2 the stau 
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Figure 6: Fitted values of tan/3 (left) and constraint from relic density combined with the pseudo- 
scalar Higgs limit (right) in the mo, ^1/2 plane after optimizing tan f3 and Aq to fulfill the 
relic density and EWSB constraints at every point. The relic density requires tan/3 « 50 
in most of the parameter space, where pseudo-scalar Higgs exchange dominates. In the 
(non-red) edges where tan (3 is lower, co-annihilation contributes. The data below the solid 
line in the right panel is excluded at 95 % confidence level. 




Figure 7: Examples of diagrams for elastic neutralino-nucleon scattering. 



becomes the lightest supersymmetric particle (LSP). However, a charged LSP cannot be the DM. As 
will be shown later in Fig. 10, for even larger values of m^ 2 this non-allowed region becomes broader, 
while the relic density requires large tan /3 away from the co-annihilation region, in which case the 
mixing in the stau sector increases and the stau becomes again the LSP. 



4 Excluded region by direct DM searches 

Scattering of LSPs on nuclei can only happen via elastic scattering, if R-parity is conserved [IJ |2]. 
The corresponding diagrams are shown in Fig. [7[ The big blob indicates that one enters a low energy 
regime, in which case the protons and neutrons inside the nucleus cannot be resolved. In this case the 
spin independent scattering becomes coherent on all nuclei and the cross section becomes proportional 
to the number of nuclei: 

- = - ? (Zf p + (A — Z)f n f (3) 

7T (m DM + m N y 

where A and Z are the atomic mass and atomic number of the target nuclei. Since the particle which 
mediates the scattering is typically much heavier than the momentum transfer, the scattering can be 
written in terms of an effective coupling. Using the notation of Ref. [53] one can write: 

^ 1 m q 27 ^-^ m q 

q=u,d,s q=c,b,t 
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where G q — AdmAq/M^ . Here M denotes the mediator, and Adm, A/ denote the mediator's couplings 

(v) 

to DM and quark. The parameters are defined by 

m pfTq = (P^qQQlP) ( 5 ) 

and similar for whilst = 1 - E g =u,d, a /t^- 

When squarks become heavy, the squark exchange diagram in Fig. [7] becomes suppressed and the 
only relevant mediator is the Higgs boson, in which case the heavier sea quarks in the sum in Eq. 
[4] become the main contributor to the elastic cross section. Remember that the coupling between 
the Higgs and down-type quarks is additionally enhanced by tan/3 [22j [M], so the charm and top 
contributions are suppressed for the large values of tan f3 in Fig. |6j 

The scalar density of the heavy quarks can be inferred from the nucleon masses, since a large density 
of heavy quarks would increase the nucleon mass. The nucleon mass can be written as Mjy = Mo+crnN, 
where Mo is the nucleon mass in the chiral limit and cr^N characterizes the effect of the finite quark 
masses. The u, d quarks contribute both as valence and sea quarks, while the heavier quarks contribute 
only as sea quarks. The strange quark content is usually parametrized by the ^/-parameter: 

_ 2(N\ss\N) 



(N\uu + dd\N) 

Phenomenologically, the ttN scattering is related to the a n N term, while the contribution of heavier 
quarks is related to the octet breaking of the baryon masses. From these data one obtains y=0.3-0.6 
[55] . Such large values are surprising, but one should keep in mind that they originate from the gluon 
splitting inside a proton, which is the non-perturbative regime of QCD and the interpretations might 
be prone to higher order corrections. 

An alternative solution for calculating the scalar quark densities in the non-perturbative regime is 
provided by lattice QCD. Such calculations allow to determine the valence and sea quark distributions 
to the g^n independently. After imposing the chiral symmetry in the lattice calculations the authors 
of Ref. [55] find that the valence quarks dominate the g^n term and the y parameter is small (< 0.05). 

The default values of the effective couplings in Eq. 5 are in micrOMEGAs [56] : = 0.033, fjfy = 
0.023, f$ = 0.26, f$ = 0.042, f!f] = 0.018, fftj = 0.26. If one takes the lower y values from 
the lattice calculations one finds [57]: = 0.020, /j?J = 0.026, f$ = 0.02, = 0.014, ffi] = 

(n) 

0.036, fj,J — 0.02. So the most important coupling to the strange quarks is reduced from 0.26 to 
0.02, which implies an order of magnitude uncertainty in the elastic neutralino-nucleon scattering cross 
section. 

Another normalization uncertainty in direct dark matter experiments arises from the uncertainty 
in the local DM density, which can take values between 0.3 and 1.3 GeV/cm 3 , as determined from the 
rotation curve of the Milky Way, see Ref. [58l [59] and references therein. 

These overall normalization uncertainties from the scalar quarks densities inside the nucleon and 
the DM densities in the Galaxy are independent of the SUSY parameter space. However, the effective 
couplings vary strongly inside the SUSY parameter space. 

The scattering cross section is proportional to the product of gaugino and Higgsino components 
squared [221 G2L so it rises rapidly, if the Higgsino component increases, which is the case for large 
values of mo, as shown before in Fig. [3J To get conservative estimates for the excluded regions, we take 
the lowest possible values of the local DM density and the couplings. The difference in couplings is the 
largest contributor to the uncertainty and the difference in excluded regions is demonstrated in Fig. 
[8] (left panel). The dotted line corresponds to the default micrOMEGAs values of the couplings based 
on the (JftN term from data, while the dashed-dotted line indicates the 95% excluded region using the 
smaller and hence more conservative couplings from lattice gauge theory. The colours indicate the 
Ax 2 = X 2 ~ Xmin contribution from the electroweak constraints: b —> X s j = (3.55 ± 0.24) • 10 -4 [60J, 
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g-2 [61J , -> y+yT < 1.1 • 10~ 8 [62J, B u ^rTv = (1.68 ± 0.31) • 1CT 4 [60], m h > 114.4 GeV [63]. 
Xmin — 2.5 in this plot for two degrees of freedom, if we consider only the measured values (excluding 
limits) minus the two fitted SUSY parameters Ao and tan/3. 

If we combine all constraints one finds that m^ 2 below 400 GeV is excluded for all values of mo, 
as shown in the right panel of Fig. [8j At low values of mo the direct searches dominate the limit, as 
shown previously in Fig. [4j At intermediate values of mo the relic density constraint dominates the 
limit, as shown in Fig. [6j while at large values the limit is dominated by the XENON100 data (Fig. 
[8)left panel). 

It should be noted that these combined limits from the LHC, WMAP and XENON100 are stronger 
than the limits from the LEP Higgs limits, electroweak and b-physics observables (red area in Fig. 
|8| left), so these hardly play a role anymore. Only g-2 yields a mild, but insignificant preference 
for low values of mo, especially if one combines the non-Gaussian theoretical errors linearly with the 
experimental errors. This will be discussed in more detail below. 

4.1 Discussion on g-2 

The theoretical value of g-2 has been reviewed in Ref. [64J, which is still in agreement with the latest 
values from Ref. [65]. We use the difference Aa M in g-2 between the experimental and theoretical value 
as input to the x 2 value and use Aa M = 302 ± 63(exp) ± 61 (theo.) • 10 -11 , so the theoretical and 
experimental errors are similar. The theoretical error from the light-by-light scattering contribution 
alone is 26, so this error dominates, but it has certainly not a Gaussian distribution. For non-Gaussian 
errors, especially errors with an equal probability in a certain interval, it is more conservative to add 
the errors linearly, as can be tested by simply comparing the convolution of two Gaussians with a 
Gaussian and a "flat top" Gaussian, where the flat region represents the interval with a constant 
probability. In the first case one finds the usual result that the probability distribution from the 
combined uncertainties equals a Gaussian with the errors added in quadrature, while in the second 
case the probability distribution equals a Gaussian with an error closer to the linear addition of the 
individual errors. To be conservative we have added the theoretical and experimental errors linearly, 
which increases the error for g-2 from 88 to 124 in units of 10 -11 , so the x 2 contribution is reduced by 
about a factor two, but there is still a preferred region with lower x 2 5 as shown in Fig. [9J left panel. 
To be consistent we have added theoretical and experimental errors linearly for all other observables 
as well. 

If we exclude g-2, the remaining observables are insensitive to the region with large SUSY masses, 
as can be seen from the flat x 2 distribution in the right panel of Fig. [9J which can be compared with 
the right panel of Fig. |HJ where g-2 was included. One observes that the XENON100 limit excludes 
m i/2 < 300 GeV for large mo, but if combined with g-2 the exclusion goes up to 400 GeV in this 
region. 

5 Summary 

If one combines the excluded region from the direct searches at the LHC (Fig. [4j left) and the relic 
density from WMAP combined with the already stringent limits on the pseudo-scalar Higgs mass 
(Fig. [6]) with the XENON100 limits (Fig. [8j left) one obtains the excluded region shown in the right 
panel of Fig. |8j Here the g-2 limit is included with the conservative linear addition of theoretical 
and experimental errors. One observes that the combination excludes m^ 2 below 400 GeV for all 
values of mo, which implies the LSP has to have masses above 160 GeV within the CMSSM from the 
constraints considered. The gluino mass has to be above 1 TeV in this case. 

As discussed in Sect. [2| the LHC becomes rather insensitive to the large mo region because of the 
decreasing cross section for the production of strongly interaction particles and the large background 
for the production of gauginos. However, in this region one obtains an increased sensitivity above 
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Figure 8: Left: A\ 2 — X 2 ~ X m in distribution in the mo,?^^ plane after imposing the electroweak 
constraints in comparison with the XENON100 limits [19] on the direct WIMP-nucleon 
cross section for two values of the form factors (dotted line: tvN scattering, dashed dotted 
line: lattice gauge theories). Values with A% 2 = X 2 ~ Xmin > 5-99 are excluded (red 
region). x^ n =2.5 for two degrees of freedom. Right: as left after imposing the combined 
constraints from LHC, WMAP, electroweak constraints and XENON100. X m in=3.6 for two 
degrees of freedom, which is slightly higher than in the left panel because of the increased 
contribution from g-2, see Fig. |9| 



the LHC limits from relic density and direct DM searches for two reasons. First, the relic density is 
inversely proportional to the annihilation rate, which leads to an annihilation cross section of the order 
of a few pb. Such a large cross section can only be obtained by annihilation via pseudo-scalar Higgs 
boson exchange in the s-channel in most of the parameter space. The mass of the pseudo-scalar Higgs 
boson can be tuned everywhere in the mcm^ plane by choosing the correct value of tan/3, which 
has to be around 50 in this case. Such a large value of tan (3 leads to a large pseudo-scalar Higgs 
boson cross section at the LHC, since it is proportional to tan/3 2 . The LHC limit on the pseudo-scalar 
Higgs mass leads to a lower limit on m^ 2 of 400 GeV for intermediate values of mo, as shown in Fig. 

El 

Secondly, the cross section for direct scattering of WIMPS on nuclei has an upper limit of about 
10 -8 pb, i.e. many orders of magnitude below the annihilation cross section. These cross sections are 
related to each other by the diagrams and kinematics. The many orders of magnitude are naturally 
explained in Supersymmetry by the fact that both cross sections are dominated by Higgs exchange 
and the fact that the Yukawa couplings to the valence quarks in the proton or neutron are negligible. 
Most of the scattering cross section comes from the heavier sea-quarks. However, the density of these 
virtual quarks inside the nuclei is small, which is one of the reasons for the small elastic scattering 
cross section. In addition, the momentum transfer is small, so the propagator leads to a cross section 
inversely proportional to the fourth power of the Higgs mass. At large values of mo EWSB forces the 
Higgsino component of the WIMP to increase and consequently the exchange via the Higgs, which has 
an amplitude proportional to the bino-Higgsino mixing, starts to increase. This leads to an increase 
in the excluded region at large tuq. 

As mentioned in the introduction, several groups have performed similar analysis. Our results are 
closest to the one of Ref. [31J. They define the limits within the frequentist approach in a similar 
way as we do. Instead of optimizing the SUSY parameters by calculating the observables during 
execution of the program, they prepare a large data base of randomly sampled SUSY points with 
the observables calculated at each point. The results for small values of m^ 2 are similar to ours, 
although they use cross sections for the XENON100 results with couplings intermediate between the 
tvN scattering and lattice gauge theories, while we used conservatively the lattice gauge theories. The 
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Figure 9: Left: A\ 2 — X 2 ~ X m in distribution of the g-2 observable alone under the constraint that 
tan /3 and Aq are still fixed by all other constraints. One observes a shallow increase of the 
X 2 value for large SUSY masses, because g-2 prefers light SUSY particles. Right: the total 
Ax 2 distribution without g-2 constraint. One observes that all points above the excluded 
region (solid line) are equally probable. Note that the combined limit is slightly reduced 
at large values of mo in comparison with Fig. [HJ right panel, while g-2 still contributes, 
even if the errors are added linearly. 



difference between tvN scattering and lattice gauge theories has been displayed in the left panel of 
Fig. [8j They display results up to m^ 2 = 2500 GeV, since they find excluded regions above this 
value, which is due to the relic density constraint [66J. In our case we do find good solutions and no 



excluded region is found above m^ 2 = 400 GeV, as shown in Fig. 10, left panel. This is probably 



due to the fact that in this region tan /3 and Aq are highly correlated, so they can be easily missed in 
randomly chosen SUSY samples. The strong correlation is shown in the right panel of Fig. 10 and 
the best solutions are obtained close to the white stripes at the top and bottom, which are near the 
stau co-annihilation region. In the white region the stau is the LSP. As shown in the right panel of 
Fig. 10 there is no preferred region above m^ 2 = 400 GeV, if g-2 is excluded and the region where 
the stau becomes the LSP is ignored. The preferred minimum for g-2 (around mo = 400, m^ = 200 
GeV (Fig. [9] left) is already excluded by the LHC data and the slight preference above m^ 2 = 400 
is solely due to the shallow tail in the x 2 distribution of g-2 (Fig. [9J left panel). How strong this 
preference is depends then on the treatment of the errors of g-2. As argued above the theoretical 
errors of the light-by-light scattering dominate and are certainly non-Gaussian, in which case a linear 
addition of the experimental and theoretical errors is the more conservative approach, so we do not 
think the preference by the region selected by g-2 and the corresponding preference for the expected 
SUSY masses is worth emphasizing in contrast to Ref. |31j . 

Our results differ significantly from results using Markov Chain Monte Carlo sampling. E.g. in 
Ref. [32] values for intermediate values of mo are excluded, which is the region of large tan /3 (see 



Fig. |6j left panel). Here the parameters tan/3 and Aq are highly correlated again (Fig. 10, right 
panel) and finding the correct minimum depends strongly on the stepping algorithm, e.g. stepping in 
the logarithm of a parameter is different from stepping in the parameter ("prior dependence"). Such 
dependence on sampling techniques largely disappears in our multistep fitting technique, since for 
each point of the mo,m^ 2 grid a unique solution is found independent of the minimzer used, so the 
frequentist approach with x 2 minimization yields the same results as a likelihood optimization with a 
Markov Chain sampling technique. 

If one combines the limits from the direct searches at the LHC, heavy flavour constraints, WMAP 
and XENON 100 using the most conservative assumptions of linear addition of theoretical and exper- 
imental errors and the lowest local relic density and matrix elements for the XENON100 limit we 
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Figure 10: Left: A\ 2 distribution of all constraints up to m^ 2 = 2500 GeV, showing that the x 2 
does not increase because of the relic density for large neutralino masses in contrast 
to Ref. [31J. The white region in the top left corner is excluded because the stau is 
always the LSP. The red region in this corner is excluded by the relic density constraint 
requiring large tan/3, which in turn cause a large mixing in the stau sector leading to 
the stau becoming the LSP again. Right: tan/3 and Ao are strongly correlated for large 
neutralino masses by the relic density constraint; here mo = 2500 GeV and m^ 2 = 2000 
GeV has been chosen as an example. 

exclude values of m^ 2 below 400 GeV in the CMSSM independent of mo, which implies a lower limit 
on the WIMP mass of 160 GeV and a gluino mass above 1 TeV. This is the CMSSM parameter region, 
where SUSY can be found. 

Note that the CMSSM relates electroweak gauginos and gluinos by requiring gaugino mass unifi- 
cation. However, the squark and gluino masses are related by the squark self-energy diagrams, which 
are valid in any model with couplings between squarks and gluinos. So the gluino limits from Fig. [5] 
of 620 GeV are model-independent. 
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